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Two aspects of isospin breaking in the decay —>■ 7r°7r°e^ W are studied and discussed. The 

first addresses the possible influence of the phenomenological description of the unitarity cusp on 
the extraction of the normalization of the form factor from data. Using the scalar form factor 
of the pion as a theoretical laboratory, we find that this determination is robust under variations 
of the phenomenological parameterizations of the form factor. The second aspect concerns the 
issue of radiative corrections. We compute the radiative corrections to the total decay rate for 

—>■ 7r°7r'^e* W in a setting that allows comparison with the way radiative corrections were 

handled in the channel —>■ 7r’''7r“e* W. We find that once radiative corrections are included, 

the normalizations of the form factor as determined experimentally from data in the two decay 
channels come to a better agreement. The remaining discrepancy can easily be accounted for by 
other isospin-breaking corrections, mainly those due to the difference between the masses of the up 
and down quarks. 


I. INTRODUCTION 

The program of analysing Ki 4 decays of the charged kaon conducted by the NA48/2 collaboration at 
the CERN SPS has so far been very successful. In the tt+tt” channel of the electron mode, £ = e [the decay 

—>■ 7r“'"7r zzj will henceforth be refered to as ], it has led, besides a more precise determination of 

the corresponding branching ratio and hadronic form factors [l[, to a very accurate determination of the tttt 
5'-wave scattering lengths Oq and ii, that constitutes a stringent test of the QCD prediction obtained 
within the framework of chiral perturbation theory M- 

More recently, the results concerning an analysis of the data obtained in the channel of the electron 

mode [the decay —>■ 7r°7r°e^ W will henceforth be referred to as AT™] have also become available Q. 
Although the number of events is lower [^ 6.5 • 10“* events in the AT™ mode vs. ^ 10® events for this 

allows for some cross checks at the level of the structure of one of the form factors, that is identical for the two 
channels in the isospin limit. The normalisation of this common form factor, as measured in the two channels, 
reads [HQ 


\Vus\fs[K+-] = 1.285 ±0.001stat±0.004y.t±0.005ext, 

(1 + SEM)\Vus\fs[K°2] = 1-369 ± O.OOSstat ± 0.006sy.t ± 0.009ext. (I.l) 

Ignoring for the time being the correction factor Sem (we will discuss radiative corrections below), the difference 
of the two values, as compared to the value measured in the channel, amounts to 6.5% in relative terms. 
This might be considered as a small difference, but given the uncertainties, it is, in statistical terms, quite 
significant. Adding all errors in quadrature this gives 


(1 -I- 6em) 


fs[K°J] 


1.065 ±0.010. 


( 1 . 2 ) 
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It seems difficult to ascribe a variation of 6.5% to the radiative correction factor Sem alone. While in 
some regions of phase space radiative corrections can reach the ±10% level, they usually sum up to ±1% in 
the decay rate. The radiative corrections to the decay mode have been discussed in several places [M 
at one-loop precision in the low-energy expansion. But no comparable study has been done for the it'™ decay 
mode. There exists an older, less systematic, analysis [l^ that covers the corrections due to virtual photon 
exchanges and real photon emission, which could provide the relevant contributions at a first stage, but its 
practical use is somewhat limited, since the expressions given there are not always very explicit, and moreover 
need to be checked. Furthermore, not all radiative corrections occurring in the charged channel 
have been taken into account in the analysis of the experimental data. These additional radiative corrections 
could affect fs[K°2] and fs[K^f] in different ways, and make up for another part of the discrepancy. 

If Sem alone does not explain the discrepancy (11.21) . one has to look for other sources of isospin-breaking 
effects. These can be due to the difference between the up and down quark masses mu and conveniently 
described by the parameter R, with 1/R = (rud — mu)/{mg — mud), where mg is the mass of the strange quark, 
whereas mud denotes the average mass of the up and down quarks, mud = {"mu ±md)/2. For instance, at lowest 
order in the chiral expansion, one has iEl 


fs^ / 3 \ 

fg[Kt,-] \ 2 r)- 


(1.3) 


Barring contributions of higher-order corrections, values of R as small as 1^ R = 35.8(1.9)(1.8) can account 
for about two thirds of the effect in Eq. (iLlt. 

Finally, there are also isospin breaking effects induced by the mass difference between charged and neutral 
pions. Most notable from this point of view is the presence of a unitarity cusp Q in the form factor describing 
the amplitude of the iF™ mode. The interpretation of this cusp is by now well understood, and as in the case of 
the 7r°7r°7r^ decay [IMl, it arises from the contribution of a 7r“'"7r intermediate state in the unitarity 

sum [for a general discussion of the properties of the Kei form factors from the point of view of analyticity and 
unitarity, see Ref. and references therein]. 

This cusp contains information on the combination Oq — Ug that describes the amplitude for the process 
7r+7r“ —>■ at threshold. Although this information probably cannot be extracted from the A'™ data in 

a way as statistically significant as the determination from the cusp in —>■ [^, it is nevertheless 

important to include a correct description of this cusp in the parameterisation of the form factor used to analyse 
the data. This necessity has been demonstrated in full details in the case of the —?► TT^n^n^ decay, and it 
is to be expected that the same attention to these matters should be paid also in the analysis of the AT™ data. 
Failure to do so may introduce a systematic bias which would make the comparison with the information on 
the form factor extracted from the AT^” data spurious to some extent. 

It is the purpose of the present note to address some of these issues. In a hrst step, we investigate the 
possible influence that various parameterisations of the form factors could have on the outcome of the analysis. 
In order to control inputs and ouputs fully, we choose to work with a simplified model, where the exact form 
factors are known from a theoretical point of view, and where one can assess the effects of various choices of 
parameterisations for the form factors used in order to analyse the numerically generated data [that we will 
henceforth refer to as pseudo-data]. This framework is provided by the scalar form factors of the pions, defined 
as 


{Tr°{pi)n°{p 2 )\m{uu + dd){0)\n) = +Ff{s) [s = {pi + p 2 )‘^], 

{Tr+{p+)Tr-{p_)\m(uu + dd){0)\n) = -F^{s) [s = (p+± p.)^]. (1.4) 

Expressions of these form factors, with isospin-breaking contributions due to the difference of masses between 
charged and neutral pions included, have been recently obtained in [2^ up to and including two loops in the 
low-energy expansion. We will use these expressions in order to generate pseudo-data, which we can then submit 
to analysis, using various parameterisations for the form factors, inspired by those in use for the analyses of the 
Ar^“ and AT™ experimental data. The reason for working with the scalar form factors is at least twofold. First, 
the form factors, with isospin-breaking effects included, are known at two loops in both channels, whereas in 
the A'e 4 case, only the form factors in the channel with two charged pions have been studied at the same level 
of accuracy as far as isospin-breaking corrections are concerned [2^ [see Ref. [2^ for a systematic study at one 
loop]. Second, the A'e 4 form factors depend on two more kinematical variables, besides the di-pion invariant 
mass. The scalar form factors depend only on the latter, and offer therefore a simple kinematical environment, 
so that the issues we wish to focus on can be addressed without unnecessary additional complications. 

In a second step, we address the issue of radiative corrections to the total decay rate of the decay 
—>■ W. Our intent here is not to develop a full one-loop calculation, at the same level of precision 

as those that exist for the decay channel into two charged pions [2-[ll|. We rather aim at providing a simple 
estimate for the radiative corrections to the total decay rate, much in the spirit of Refs. [l2| and |24l [ or, on 
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a more general level, of Ref. [2^. This will allow us to assess how much of the discrepancy m has to be 
ascribed to other isospin-breaking effects in the form factors, such as discussed above. 

The remainder of this study is then organised in the following way. First, we give (Sectionln]) a theoretical 
discussion of the structure of the scalar form factors of the pions using the explicit expressions obtained in Ref. 
[ 2 ^. We will thus adapt the discussion of Ref. Q to the case at hand. Working on this analogy will allow 
us to give an assessment of some additional assumptions regarding the structure of the form factors implicitly 
made in Ref. Q. Next, we generate pseudo-data ISection llllll using the known two-loop expressions of the form 
factors, that we then analyze using various phenomenological parameterisations, that do not necessarily comply 
with the outcome of Section El The purpose here is to discuss in a quantitative way the possible systematic 
biases that can be induced by these different choices. The last part of Section Hill addresses the determination 
of the combination Ug — Ug of S'-wave scattering lengths. Radiative corrections, aiming at an estimate of the 
correction factor Sem in Eq. (EB, are discussed in Section HVl We first compute radiative corrections to the 
decay rate in a similar set-up to the one used for the treatment of the data in the channel, in order to 
obtain a meaningful comparison between the two channels. Then we compute the effects of additional photonic 
corrections, not included in this treatment. Finally, we end our study with a summary and conclusions. Two 
Appendices contain technical details relevant for the discussions in Section El and in Section llYl respectively. 


II. DESCRIBING THE CUSP: THEORY 


According to the general analysis of Ref. [13 , the occurence of both and tt+tt intermediate states 
at different thresholds leads to the following structure for the scalar form factor of the neutral pion Fg (s) [M-j^ 
stands for the charged-pion mass, whereas Mt^o is the mass of the neutral pion]: 


e 




•> 

Ff (s) + Ff (s) < s < 4M2] 


(II.l) 


where F^ (s) is a function of s that is smooth as long as no other threshold, corresponding to higher intermediate 
states, is reached. Here 5{s) represents a phase. It can be chosen arbitrarily, as long as it is also a smooth 
function of s. The cusp at s = observed in the differential decay rates corresponding to this simplified [as 
compared to K^a] situation then results from this decomposition, since 


F^s (s) 


2 



-^o”(s) 

2 

+ 

2 

Ff{s) 

2 

— 2Im 

2 

V(s)-^r"*(s)' 




Ff{s) 

-f 2Re 

Ff{s)Ff*{s) 


[s > 4M2] 


[dM^o < s < 4M2] 


(II.2) 


Apart from the dependence with respect to the second kinematical variable Se, the empirical parameterisation 
used for the fit to the AT™ data, Eq. (9.1) in Ref. Q, complies with this general representation provided [the 
variable S-,^ used in this reference corresponds to the variable s used here]: 

1) Fi (s) is parameterised as a real constant times o'(s), with 


d-(s) = 



/ 


l + q^ 


(11.3) 


for 4M^o < s < 4A1^, or < 0, with s = 4M^(1 -|- q^). 

2) Fl° (s) is set to zero (its value for s = 4M^) for s > 4M^ (q^ > 0) 

3) For 4M^o < s < 4M^ {q^ < 0), Fq° (s) is replaced by a constant, equal to Fq° {4M^). 

A more theoretically based parameterization, adapted from the simple discussion of the cusp in —>■ 

given in Ref. [l^, is considered in Sec. 9.4 of Ref. Q, though not used for the data analysis. As compared 
to Eq. (III. 21) . its validity also rests on additional assumptions, which, once transposed to the present situation, 
read: 


1’) The phase 5{s) can be chosen such as to make the two functions Fq (s) and Fi (s) simultaneously 
real, so that Eq. (III. 21) takes the simpler form 






[s > 4M2] 

[4M^o <s< 4M2] 


(II.4) 
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2’) is related to the scalar form factor Fg{s) of the charged pion, multiplied by a combination of 

the two S'-wave tttt scattering lengths Oq and Og in the isospin limit, 

= (n.5) 

In view of the discussion in Ref. Q, Fg{s) should be identified with the phase-removed scalar form factor 
of the charged pion. The latter is given by Fg {s), where the phase ^J(s) is defined as Fg{s -I- fO) = 

e2*'5o(®)F^(s-fO). 

Our purpose in this Section is twofold. First, we will rewrite the two-loop representation of the form 
factor Fg (s) obtained in Ref. [2^ in the form (III. 21) . that makes the cusp structure explicit. Second, we will 
assess to which extent the additional features mentioned above and assumed in Ref. Q are actually reproduced 
by the structure of the form factors at two loops in the low-energy expansion. In particular, we will establish 
the precise relation between Fg{s) and Fg{s) in Eq. (III. 51) at this order. In what follows, and unless otherwise 
stated, it will always be understood that s > Furthermore, in practice s > will actually mean 

dMj < s < M|-, where Mk is the mass of the charged kaon, so that we need not worry about thresholds other 
than those produced by two-pion intermediate states. 


A. The cusp in the one-loop form factor 

We start with the study of the cusp using the one-loop expression of the form factor Fg°(s), 


(s) = F^ (0) 


I + Oo s + IGtt- 


-Jo{s) 


- 167rF^(0)v3g(s) J(s). 


(II.6) 


In this expression, a'g denotes a subtraction constant, that we need not specify further for the time being. The 
loop functions Jo(s) and J(s) are given by 


>/o(s) = 

J» = 


s 


f°° dx I 
/- rxMx) 

74^2 X X — S — lU 


r°° ^_ i_ 

^ X — s — iO 


7Tcr(a;), 


with 


. X / , X / 4M2 

cro(s) = \ l -, cr(s) = \1 -. 


(II.7) 


(II.8) 


The functions ‘/?o°(s) and iPq{s) denote the lowest-order real parts of the S-wave projections of the amplitudes 
of the processes 7r°7r° —>■ 7r°7r° and 7r°7r° —>■ tt+tt”, respectively. Their expressions read 


00 X X 2 , X X , s — 4M2 

— ®00; x'®) — T i 


(II.9) 


with 


aoo = ^ («o + 2ao) (^1 “ (aj] - ag) -b (2aSS - 5ag) (II.IO) 

and = M2 - 

In the range of s under consideration, the function Jo(s) is complex, but both its real and imaginary 
parts are smooth. 


Jo{s) = [2 + cro(s)Lo(s) -b i7rcro(s)], Lo(s) = In 


1 - go(s) 
1 + cro(s) 


[s > dM^o], 


(II.ll) 


whereas J(s) may be rewritten as 

J(s) = J[0](s) + J[l](s) X 


—ia(s) [s > 4M2] 

+d(s) [4M2 o < s < 4M2] 


( 11 . 12 ) 
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The two functions and j[^l(s) are smooth, and read 




IGtt^ 


2 + a{s)L{s) , jW(s) = - —, 


where 


L{s) = In 


1 - cr(s) 


[s > AM^o]. 


+ cr(s). 

In these expressions, the definition of cr(s) has been extended below s = 4M^ byHi 


s) = < 


4M^ 

1 -^ = <t(s) [s > 4 M 2 ] 


AM£ 


— 1 = ia{s) [AM^o < s < AM^] 


(11.13) 

(11.14) 


(11.15) 


According to Eqs. (III. Ill and (111.21) . J{s) exhibits a cusp structure at s = AM^. One thus obtains, at this order, 
the decomposition of the form (III. II) for Fg (s), with 


^0 (■s) — -P's (0) % 1 + oj s + 


0. , 


2tt 


[2 + ao{s)Lois)]}~F^iO) 


V>o{s) 


2 + a{s)L{s) +0(^6), 


is) = Fs-(0)^g(s)d(s) + 0{E^), 


(11.16) 


provided one factorises the global phase 

^(s) = ^^o(s)‘^r(s) (11.17) 

Therefore, up to so far unspecified higher order corrections, the one-loop expression of the form factor can be 
brought into the form (III. 111 . Both functions Fq (s) and F^ (s)/(t(s) are real and smooth for s > AM^o at 
this stage. However, the expression for J^f°(s) in (III.161) does not quite comply with Eq. (III. 51) . Whereas at 
this stage Fg{s) is equal to the constant Fg{0), which, at this order, can be identified with the phase-removed 
form factor, the combination of scattering lengths that occurs in Eq. (111.51) corresponds to ‘Po{AM^), where 
7’o(s) is the expression of </5g(s) in the isospin limit. Thus, at this order, the expression (III. 51) misses both the 
dependence with respect to s in ‘Pq(s), and the isospin-breaking corrections in the scattering lengths. 

For later reference, we briefly extend the discussion to the scalar form factor of the charged pion. At one 
loop, it is given by 


F^(s)=F^(0) 


Zg s + 167r(pg (s)J(s) 


- 16^E/(0)i^g(s)Jo(s). 


(11.18) 


Besides the subtraction constant Og, that differs from ag° (they become identical in the isospin limit), this 
expression involves the lowest-order real part of the S-wave projection of the amplitude for the scattering 
process 7r“''7r“ —>■ 7r'''7r“, 


s — AM^ 

Po (s) = a+- + K- ’ 

(11.19) 

where 


«+- = 2 (2^0 + a-o) “ ^+- ^ ^ 

(11.20) 

After having factorised the global phase 


5{s) = -\ao{s)^^,{s)^^+0{E% 

(11.21) 


one can decompose Fg (s) according to Eq. (III.1|) , with 


1 + + 


‘Po (s) 


2 -I- a{s)L{s) [■ — Eg (0) 


27T 


[2 + ao(s)Lo(s)] + OiE^), 


= -Eg-(0)(^+-(s)d(s) +0(E6). 

Both functions Fq{s) and Fi{s)/a-{s) are real and smooth for s > AM'^o at this stage. 


( 11 . 22 ) 
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B. The cusp in the two-loop form factor Fg (s) of the neutral pion 


Let us now go through the same analysis, but with the two-loop expression of the form factor. The 
expressions of the pion scalar form factors at two loops and in presence of isospin breaking have been worked 
out in Ref. using a recursive construction based on general properties like relativistic invariance, unitarity, 
analyticity, and chiral counting. The scalar form factor of the neutral pion can be written as(^ 


(■ 5 ) — P's (0) (^1 + aj s -I- 65 


- 167rF^(0) (pg(s) 


^0 


1 


Mi 


+ ^ Ff{0) ^il\s)Ms) + eoris)Kns) + 2eons)K!iis) + eoris)K3is) 

+ do"\s)Kns)+dT\s)K^is) + [16^- 2] Ms) 

A/f4 _ 

-2-^FS{0) eHs)j{s)+2iM\s)KM)+e'>mns)+6^Hs)K^{s) 

7T 

+ 2^M\s) [WJo(s) - 2] J(s) + Ai^,(s)£"(s)l + 0{E^). 


1 -I- Og S -I - ifQ (s) 

TT 

1 -I- Og S -I - ipQ (s) 

TT 

tim, 


Jois) 

As) 


j(2;o)/ 


f(3;0). 


( 11 . 23 ) 


In this formula, the functions ioo^(s),.... A^^s),... are polynomials of at most second order in the variable s. 
Their expressions can be found in Ref. [2^, except for Aiia;(s), that reads 


Aii,(s) 8 ^2 ^+ 0^2 
It is also useful to be aware of the relations 


9Mi 


f Mio 

^+0 - + 26+0 f 1 + 


( 11 . 24 ) 


^^oA\s) = 2§-Ml\s)]^ £_M){s)=2^y,\s)Ms). ^'iA\s)=^%[Ms)t e^\s) = ^^^Ms)A-{s). 

-^^TT -^‘^TT -^^TT 

( 11 . 25 ) 

In order to achieve the decomposition (III. II) . we need to extend the decomposition of the function J{s) 
in Eqs. (III.121) and (III.131) to the other functions, denoted generically by K^{s), that appear in the expression 
(III.231) . This may be done as follows. First, we may observe that, like Jo(s) or J(s), these functions can also 
be defined by a dispersive representation of the form [for Jo(s) = M(s)y one has ^[[(s) = cro(s)/167r, whereas 
for J(s) = M(s), ko(s) = cr(s)/167r, see Eq. (III. 71) ] 


P:(s) = - / — 


1 


dx 


X X — s — iO 


KM- 


( 11 . 26 ) 


Explicit expressions for the functions A:“(s) are given in Appendix]^ For the set of functions K^{s) and (s), 
one has Sthr = 4M^o- These functions will therefore each develop an imaginary part for s > Sthr = 4M.n.o, 
Im.A“(s) = fc“(s)0(s — 4M^o), while the real part displays a cusp at s = 4M^o, but is smooth for s > AM^^o. 
The situation is different for the remaining functions, Kn{s), K^{s), and iC^{s), for which sthr = 4M.^, so that, 
in a generic way, they have the following structure 


KZ{s) = Rei?“(s) + 


M[s) [s > 4M2] 


0 


= Rei?“(s) + :^x 
cr(s) 


[AMlo < s < AMl] 
id{s) [s > 4M^] 

0 [AMI, < s < 4M2] 


( 11 . 27 ) 


In general, the function k“{s)/(T{s), although real, is not smooth for the whole range s > AM^,, but only for 
s > AM^. Suppose one can find a function fc“(s) such that it coincides with fc“(s) for s > AM'^, and such that 
k^{s)/a{s) is real and smooth for all s > AM"^,. Then one can perform the decomposition 


K{s) = xM(s)-kXW(s) X 


—id{s) [s > 4M^] 

+d(s) [AMl, < s < AMl] 


( 11 . 28 ) 
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in terms of two real and smooth functions i4r[°l(s) and given by 




ReK“(s) 


[s > 4M2] 


ReK^{s)-ik° (s) [4M^o < s < AIVQ 


= - 


Kjs) 

a{s) ■ 


(11.29) 


Such a decomposition can indeed be achieved for the various functions considered here, as discussed in detail 
in App. |Al The decomposition (III. Ill of the form factor now follows immediately, with 


J-J” (s) = Ff{0) (l + afs + bfs^) 


TSttFI 

-167rF^(0)^S(s) 


1 + Og s H- ipg (s) 


l + aJsH (pg (s) 


Ms) 

Mis) 


+ If + dTmM) + Mo"\s)K-is) + 

M TT 

+ - 2] Ms) 

- 2 || F^iO) j™(s) + e\s)Kf\s) + 2e’^Hs)M°\s) + e\s)Kf\s) 

M TT 

+ 2^^M\s) [16^2 Jo(s) - 2] M\s) + + OM), (11.30) 


and 


e'‘^(")-^r"(s) = -<5-(s)i 167rAJ(0)v?S(s) 


1 + 05 S H- (fg (s) 


JW(s) 


-2^Ffm^^g’^\s) X 16^2Jo(s)JW(s) 

M TZ 

+ 2 S F^iO) k(°)(s) JW(s) + MHs)kM{s) + 2eMHs)M^\s) + eHs)Kf\s) 

M TZ 

+ 2^f’^\s) [I67r2 Jo(s) - 2] .M (s) + AiUs)M^^\s)] ] + 0{E^). (11.31) 


Both functions Fq° (s) and Fi° (s) / a{s) are smooth for s > 4M^o, but complex. It remains to discuss 
the phase Sis). If we want to make the function Eg (s) real, while keeping it smooth, then its choice is unique. 


Hs) = ^Ms) 


Mi^) + Mis) 


(11.32) 


with "00 (s) given by 


■Ms)% is) = 


Ml 

F^ L 


is)kois) + ^oo°\s)kiis) + dM\s)ki is) + ^m°\s)k2is) 


Hoo°\s)k3is) + Mo"'\s)kMs) + Mo^\s) ^cro(s)cr(s)L(s) 


Stt 


(11.33) 


Note that ^q°(s) differs from the quantity t/'o°(s) defined in Eq. (4.6) of Ref. [l^ by the presence of the function 
Lis) instead of L(s) in the last term between square brackets, see Eq. (IA.2E This makes tTo(s)^o°('®) ^ smooth 
function for s > 4 M^o, whereas ct(s)T(s), and hence tpg’^is), displays a cusp at s = 4 M^. Making use of Eq. 
(III.271) . the removal of the phase Sis) indeed leads to a real and smooth expression for the function J'g”(s): 


Fq is) — Fg (0) (1 + Og s + 65 


+ 87TFfiO)<p^g°is) 

-m7rFM0)Ms) 


1 + Oo s 


Re Jo(s) 


1 + 05 S H- (fig (s) 

TT 


jM(s) 































Fi 


+ ^ FfW te(s)Re Jo(s) + + 


■^00 

>( 2 ; 0 ). 


+ e;ir(s)ReF2'(s) + ^(s)ReF3^(s) + KJ i^) 


(3;v), 


Mi 


- 2 F^{0) jm(s) + + 2^i^-^\s)K^°\s) 

TT ^ 

+ ^i^'>{s)Kf\s) + A,Us)iC^^°Hs) 

+ 8Ff{0) [(^g(s)]2 (l67r2jM(s) - 2 ) Re Jo(s) 
-8F^(0)(^S(s)(^™(s) [I6^2j^e Jo(s) - 2] J[°l(s) 


-Ff{0) K(.)]^ 


2ReF^(s) + 8Re Jo(s) + - 1 - 


4Mio 


+ C>(F® 


(11.34) 


As far as Fi (s) is concerned, we may even proceed in a more direct way by noticing that, up to higher order 
corrections, Eq. (IIL31I) rewrites as 

xFs-(0)|l + a5s + 16^<^+-(s)J[°l(s)-87r:^^^S(s)Jo(s)| + 0{E^), (11.35) 


with [for the notation, see Appendix [A] 


rois) = 2 


Mi 


F^ cr(s) 


6°Hs)ko{s) + ^^^’^\s)k 2 {s) + ^i^\s)k^{s) + {s)k^{s) 

+ ^i^’°\s)k^{s) +Ai^^{s)k^{s)\. 


(11.36) 


Now, the phase that appears factored out on the right-hand side of this equation can be identified with the 
phase S{s) on the left-hand side, since the difference generates contributions of order 0{E^), that are neglected 
anyway. Taking into account Eqs. (III.211) and (III.22|) . one finally obtains 


J 7 (s)=d(s) F.(s) + 0 (F«), 


(11.37) 


with 

F^{s) = e*^(")Fo"(s) + 0(F®). (11.38) 

It is possible to give a more precise interpretation of the combination </5g(s) -I- V'o (s) that occurs in (III.37p . To 
this end, let us recall from Ref. [22| that the £ = 0 partial-wave projection /g (s) for the scattering amplitude 
of the process 7r°7r° —)► tt+tt” is given, at order one loop and for s > 4M^o, by 


/o(s) = </^§W + ^oW + *‘^5(5) 




)^\ 


™(s) + a(s)^rW0(s-4M2) 


0{E^), 


(11.39) 


where '0o(s) is defined in Eq. (4.15) of Ref. [l^ [the contribution A2'4>q{s), of second order in isospin breaking, 
is numerically quite small, and is omitted for simplicity]. It differs from '!/'o(s) by the replacement of fc 2 (s) by 
^ 2 ( 5 ) in Eq. (II1.36E Then applying the decomposition (III.281) to /g (s), one finds 


( —ia{s) [s > 4M^] 

e-*'(*Vo(s) = /o"™(s)+/o™(5)x , (11.40) 

[ +(t(s) [ 4 M 2 o < s < 4 M 2 ] 


with f§^°\s) = (p^{s) + 0{E^), and = -ip^{s)ip+ {s) + 0{E^). 

We may summarise this theoretical study of the cusp in the scalar form factor of the neutral pion with 
a couple of remarks: 

• It is, in general, not possible to chose the phase (5(s) in Eq. (III.II) such as to make both Fg”(s) and Ff°(s) 
real simultaneously. A relative phase remains, see Eqs. (III. 371) and (III.38E At lowest order, this phase 
is given by the S-wave projection of the inelastic rescattering of a pair of neutral pions through a pair of 
charged pions, cf. Eq. (III.21E 
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• The structure of (s) is more complicated than just the product of the scattering length corresponding 
to this rescattering amplitude times the phase removed scalar form factor of the charged pion. At the order 
we have been working, it involves the decomposition (III.401) of the S-wave projection of this amplitude 
times the part of the decomposition (III.ll) of Fg{s). This is different from the phase-removed form 

factor, as already seen at order one loop: 


r 0 [s>4M2] 

- F^{s) = (s) X \ + 0{E^). (11.41) 

{ 1 [4M2„ < a < 4M2] 

Note however that this difference only concerns the region 4M^o < s < which contributes very little 

to the total decay rate as defined by Eqs. (IIII.ll) and (IIII.2I) below. 


C. Description of the two-loop form factor Fg{s) of the charged pion 


We now briefly address the scalar form factor of the charged pion. The issue here is not to describe the 
cusp, that occurs below the physical threshold at s = 4M^, but to provide the expressions that will be used in 
the sequel. Again, we will rely on the results obtained in Ref. [22| , and rewrite the form factor Fg (s) at two 
loops in a way that is adapted to our purposes. In particular, we will consider the phase-removed form factor, 
which reads, in the relevant domain s > 4M^, 


= Fs-(0)(l + a5s + 55s2) 


-87rF|(0)^g(s) 

-f 167rAJ(0) 

Mt p. 

r Q 


1 -f O 5 S -I-(/3q (s) 


Re Jo{s) 


1 -I- aj s -I (fiQ (s) 


Re J(s) 


^ s“(0){4°Hs)Re Ms) + (s) + 2e(2;°)(s)ReK°(s) + eKs)ReMo{s) 

TT 

+ (s)Re (s) + (s) _ 2 ] Re Jq (s) } 


Mt 


+ 2 ^ F^M i ^i°_\s(s)Re Jis) + er:ri(s)Reifi(s) + 


(i;±)/ 


j(2;±). 


Re 1 ^ 2 ( 5 ) + xt: (1 — 


32 


4Mi 


+ ^fitk^)ReKM) + (s) + Kiis)+ 




(3;a). 


+ <-;ks) 


(iGTr^Re Jo(s) - 2 ) Re J(s) 


1 


1 - 


4M2 


0{E^). (11.42) 


The functions Kt){s) that appear in this expression have again a dispersive representation of the form displayed 
in Eq. (III.261) . The absorptive parts are in part given in Appendix]^ For the remaining one, we have 


kZ°{s) = ^K{s), r«(s) = ^r(s). 


r(s) 


r(s) 


(11.43) 


for the functions corresponding to Sthr = 4M^o, and, for those whose dispersive integrals start at Sthr = 4M.^, 


ktis) = 


cr(s) 


■ Lq{s — 4A.n.) , 


k^(s) = 


3 Mi 


■Llis-AA^)., 


(11.44) 


StT (To(5 — 4A7r) ^ IGtt sa(s) ^ 

with the definitions of the functions cr(s), cro(s) and Lo{s) for s > given in Eqs. (III. 81) . (III.Ill) , and (III.15E 


III. GENERATION AND ANALYSIS OF THE PSEUDO-DATA 


In order to study the effect that particular choices of phenomenological parameterisations of the form 
factors can have on the output, we will first generate numerical data sets for the scalar form factors of the neutral 
and charged pions. The pseudo-data in question consist of the (unnormalized) decay distribution defined by 


dsdse 



4MA 


-\FS; 


data 


(s)PA3/2(M2 




s,se), 


(nil) 
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with X{x,y, z) = +y^ + z^ — 2xy — 2xz — 2yz. The total decay rate is obtained by integrating the distributions 

(IIII.ll) . convoluted with the K ^4 phase space, over the whole physical range [we now consider the electron mode 
only, and set rrie = 0]: 


rMi 


=Af 


ds 


n(M^ — 

dSf, 




(s,Se) 

dsdSf. 


(IIL2) 


Since we consider the scalar form factor instead of form factors, the integration with respect to se, involves 
the phase space only. The overall normalisation factor has been chosen to be the one of the KeA decay. 


Gl\Vus\^ 1 
3- 2i27r5^|, |f^(0)|2’ 


(IIL3) 


up to the factor l/|Fg(0)p, introduced for convenience. 


A. Form factors and input parameters used for the generation of pseudo-data 


The form factors involved in the preceding expressions are considered as known exactly and are con¬ 
structed as follows. For we will basically use the decomposition of Eq. (III. Ill , with (s) given by 

Eq. (in. 341) . and given by Eqs. (III.37|) and (III.38|) . Unfortunately, for some of the functions involved in 

Eq. (|II.34I) . like Re or KeK!^{s), explicit analytical expressions are not known. Eor a numerical approach, 

we could use their dispersive representation, as given by Eq. (III.26p . We have however found it more convenient 
to start from expressions upon which we have full analytical control. For that purpose, one may replace the 
functions ReKI^{s) by the corresponding functions Re A'°(s), and likewise Re.K'“(s) by Re.K'„(s). We also drop 
the contribution proportional to Ai^a;(s)- In the range of s we are interested in, 4M^o ^ s < M^, the difference 
induced in (s) by these changes is numerically very small. For the scalar form factor of the neutral pion, 
the resulting expression then reads 


with 


and 


f; 


S;data 


(s) = 


•^0"data(s) + •^l"data(^) ^^0 < S < 


(III.4) 




i^data('S) — Fs (0) + aj S -f 65 


-f 87rF_5 (0)(^Q (s) 1-I-aj s ReJo(s) 


-16^F|(0)<^S(s) 

Mi 




1 + as ® d—V^o ('®) 

TT 

+ Ff{0){^^l\s)ReMs) + 

+ ?Oo’ + Coo’ ^ ('®) + ?Oo’ ^('®) 


ReK^{s)+ 

ReKlis)} 


Mi 


- 2 PsiO) + 2^(2^) (s) A™ (s) + 


+ 8Ff{0) [v?S(s)]" (l67r2 jM(s) - 2 ) Re Jo(s) 
-8F^iO)y,-{s)ip°°{s) [idTT^Re Jo(s) - 2] 


+ Ff{0) 


2ReA:^(s) -h 8Re Jo(s) + „ U “ 


4Mir 


(III.5) 


-^r;data(s) = c” 4 ™ d(s)F|(0) 


X 1 -h ajs -f 


‘Po (s) 


r{s)L{s) 


F^ (0) pUs) 

Fm 


[2 + cro(s)Lo(s)] 


Poi^) + ^ ^ + 25(2^)(s)a(s)L(s) + S^Hs) L 


(III.6) 

4^) 
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Finally, in the case of ^ 5 .data('®)> expression (III.42|) of the phase-removed form factor, replacing the 

functions K'^°{s) by K^{s), and the functions K^{s) by Kn{s), respectively. This then gives 


^S;data(s) — ^si^) + “S'® + ^S®^) 


-8^F^(0K(s) 


1 -|- aj s H—tpo (s) 


Re Jo(s) 


-f 167ri^^(0)(/j+ (®) 
Mi f , 


l + Os®4 —‘Po (®) 

TT 


Re J(s) 


Mi 


r(0) 


+ 2ef:ti(®) 


+ 2cf:;i(®) 




Re K2{s) -\- 


32 


^i°(®) + 2?i^ 

«(»; 

2] Re Jo(s)} 


®) + el-fi(®) 

j Re 

/ 4M2\ 

1 


-k 


ef:ti(®)+?l-fi(®)lReK3(®) 


(ibTT^Re Jo(s) — 2) Re J(s) H-( 1 — 

^ ' 64 ' 


4Mi 


(IIL7) 


In the sequel, we will generate pseudo-data using the expressions presented in this subsection, considered to 
provide exact descriptions of the form factors. In particular, it is understood that higher-order contributions 
are considered as vanishing. As already mentioned, in order to work within a framework where we deal with 
fully analytical expressions of the form factors, we have made some approximations as compared to the two-loop 
expressions discussed in the preceding Section. Numerically, these differences are small, but most important is 
that the approximations we have made preserve the general features of the form factors as described after Eq. 
dlLdOl) . 

For the numerical generation of the pseudo-data, we need to fix the values of the various parameters that 
occur in the expressions of the form factors. As we want to compare different methods of analysis of the Kn 
form factors, we only aim at choosing values that are representative of the expected situation in these decays, 
with some limited arbritrariness in this choice. In the following, we consider the case 

ag = 0.22 al = -0.045. (III.8) 


We fix the subtraction constants of the form factors by requesting that the charged scalar form factors has the 
typical values = 0.60 fm^ and Cg = 10 GeV“^ [l^, leading to 

al = 2.63 GeV^ bl = 2.96 GeV'‘ (III.9) 


Using Ref. [l^, one can compute the isospin-breaking shift between aj” and aj. Assuming that Cg = cl°, we 
obtain 


af = 2.60 GeV^ bf = 3.24 GeV^ (III.IO) 

leading to r^o = 0.59 fm^. For the remaining parameters, we use the same input values as in Ref. We 
normalise the scalar charged form factor to unity at s = 0 and rescale the neutral one accordingly: 


FfjO) 

FliO) 


0.99. 


(III.ll) 


As an illustration, we quote the values obtained for the total decay rate of Eq. (IIII.2I) from these inputs, for 
Vus = 0.2255: 


F’^" = 0.73 • 10"22 GeV. (III.12) 

We also show, on the left panel of Fig. [U the various contributions to the form factor obtained with our 

input values. For comparison, the right panel shows the equivalent results in the case of the parameterisation of 
the K^2 form factor A4(s) discussed in Section 9.4 of Ref. 0. It corresponds, for AIi(s), to the parameterization 
of Eqs. (IIII.13I) below with /"' = 0 and the remaning parameters fg/fs and f'J/fs fixed at the central values 
given in Table 3 of Ref. Q. For A4o(s), we take the expression of Fq° (s) given in Eq. (IIII.14I) . The parameters, 
fs, fso/fs, f'so/fs, and f'Jo/fs it involves are determined from a fit to the phase-space distribution given by 
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Ps{^,s,=Q)IPsmls,=Q) 



FIG. 1: Contributions, as functions of s (in GeV^), to Fg° {s) normalised to its value at s = 4M^: (s) (black solid), 

ReT^f (s) (blue dashed), and Im(green dotted-dashed). The left panel corresponds to our representation 
(IIII.4II . (IIII.5II . and (IIII.6II of the exact form factors with the values of the parameters discussed in the text. The right 
panel corresponds to the parameterisation used by NA48/2 in Ref. for the form factor, with the parameters 
chosen as described in the text after Eq. (IIII.12II . 


Eq. (9.1) of Ref. @ at Se = 0, and with the values of the coefficients a, b, and d taken at their central values 
as shown in Table 1 of that same reference. The overall normalisation N is fixed such that the distribution is 
equal to unity at s = 4M^, Se = 0. One can observe similar features in both plots on Fig. [T] [the absence of an 
imaginary part in the right-hand plot has already been discussed at the beginning of Sec. m, suggesting that 
our subsequent analysis, based on the scalar form factors of the pion, has also some bearing on the Kf.^ form 
factor. Note that due to our choice of phase space in Eqs. (IIII.ll) and (IIII.2I) . in practice the region of interest 
in Fig. [T] corresponds to s < 0.15 GeV. 

In summary, the form factors used in order to generate the pseudo-data are defined by the expressions 
(1111.51) . (IIII.6I) and (IIII.7I) (with vanishing higher-order corrections), together with the values (central values 
for all parameters, no error bars) of the parameters specified above. The form factors thus defined will be 
referred to as the “exact” form factors, considered to represent the “truth” to which we will fit different model 
parameterisations of the form factors, in order to obtain a quantitative determination of the possible biases 
different parameterisations can have on the output of the analysis. 


B. Phenomenological parameterisations 


In order to mimic the situation in Kf.^ decays, the pseudo-data generated with the exact scalar form factors 
of the pion will now be analyzed using approximate phenomenological parameterisations. For the analysis itself, 
we will consider a framework close to the experimental set up for the [IHl] and the 01 decay channels. 
From here on, we therefore also use = s/AM'^ — 1 in addition to s, the square of the center-of-mass energy 
of the dipion system. The region below the cusp corresponds to < 0, while positive values of describe the 
region above the cusp. 

For the charged-pion form factor, this means that we consider a parameterisation of the form 


F^is) = fs 


£f £ff £fff 

l + fq^ + 

Js Js Js 


(III.13) 


In the case of the neutral-pion form factor, we consider two parameterisations: 
• Model 1: 


FS is) = /, 


sO 


l + ^q^ + ^q^ 

JsO JsO 


Fi (s) = -2/3(a[; - al)fsa{s) 


(III.14) 


• Model 2: 


(s) = fso 


1 + ^qF 

JsO 


n, 

JsO 


FI (s) = a{s)^^{s)fs 


1 + fg^ 

Js 


£ff £fff 

Js Js 


^-^<7o(s)ipa(s) 


(III.15) 
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In order to check the influence of possible higher-order terms in the expansion, as compared to the parame- 
terisation considered in Ref. Q, we have introduced a coefficient which can be set to zero or kept as a free 
variable in the fit. 

The first parameterisation with /"' = 0 reproduces exactly the one that was considered in Sec. 9.4 in 
Ref. 0 . The second parameterisation incorporates more information gathered from the theoretical discussion in 
Section|TTl while remaining sufficiently simple. Although we have chosen not to distinguish them, the parameters 
/s, /', /", and /"' appearing in Eqs. (IIII.14I1 and (IIII.15I1 are not, a priori, identical to those occurring in the 
expression (IIII.13I) . One issue of the analysis we will present is precisely to determine to which extent e.g. 
/s in Eq. (IIII.14I) should be expected to agree with /g in Eq. (IIII.13|) . Our first task is therefore to provide 
reference values for the various parameters. This is done by performing, in Eqs. (IIII.5I) and (IIII.7I) the Taylor 
expansion around = 0, thus obtaining /go, f'sQ, fJo from the former, and /g, /', /", /"' from the latter. In the 
case of Fg(s), we neglect the small half-integer powers of q arising in the expansion, which do not contribute 
significantly in the vicinity of = 0. The resulting values are shown in the last column of Table ID These 
expansions are not supposed to provide accurate descriptions of the corresponding form factors over the whole 
physical range. 

The comparison with our various fits will illustrate how different the parameters extracted from the fit 
and those describing the real Taylor expansion can be, and will thus give information on the possible bias 
introduced by the fitting procedure. For convenience, in the following fso,fgQ... will be called “neutral” 
parameters, whereas fs, fg ■ ■ ■ are referred to as the “charged” parameters. 


C. Fitting procedures 

In order to stay close to the NA48/2 experimental set up, we will thus assume that we have measurements 
of I Eg (s)p at the 12 points corresponding to the barycenters of the experimental bins, and we assign a statistical 
uncertainty derived from the number of events collected in each bin (~ 2900 events in the first two bins, and 
^ 5900 events in all the other ones), without any correlations between the bins. As the parameterisations given 
in Eqs. (IIII.14I) and (IIII.15I) depend on the S-wave tttt scattering lengths, our will also include an uncertainty 
on these quantities in order to mock up the fact that in the real analysis these scattering lengths are determined 
from the charged form factor. Here we will use the experimental information on these quantities, namely the 
latest NA48/2 combination of AT —>■ 37r and Kg 4 results Q 

ag = 0.2210 ± 0.0047 ± 0.0040 , ag = -0.0429 ± 0.0044 ± 0.0028 , = 0.92 , (III.16) 

where we combine statistical and systematic uncertainties in quadrature. One should notice that the central 
values are close (but not identical) to the “true” values used to generate our pseudo-data. 

We consider the following methods to determine the coefficients of the above models. 

• Method A: fit of |Eg°(s)p for all points to determine all the (neutral, charged) parameters (setting 
/"' = 0 to ensure a reasonable convergence of the fit), assuming the equality of the neutral and charged 
normalisation (/g = fso) 

• Method B: fit of |Eg°(s)p for all points to determine all the (neutral, charged) parameters, setting /"' = 0 
and keeping the normalisations fs and /sq distinct 

• Method C: fit of |Eg°(s)p to determine the neutral parameters, injecting information on charged param¬ 
eters by adding to the ^ contribution corresponding to a fit of the charged form factor Eg(s) to the 
polynomial expression (IIII.13I) . effectively identifying the charged parameters in the models for Eg°(s) 
with the parameters occurring in the charged scalar form factor. 

In method C, we generate pseudo-data points for the charged-pion scalar form factor with energies corresponding 
to the barycenters given in ref. [I| , and use the relative uncertainties for Fs (combined in quadrature) quoted 
for each bin in this reference, without correlations. In agreement with ref. [l|, we add an overall 0.62% relative 
uncertainty, completely correlated between all the charged bins. The curvature of the charged form factor Eg (s) 
being more pronounced than that of the scalar form factor, a g® term must be included in the polynomial 
in order to obtain a good description of the form factor over the whole kinematic range. 

We give the resulting x^in (obtained from the best-fit values of each method). Even though each model 
provides through its fit a value of fso, one can also determine the latter by considering the branching ratio. 
In this case, fso is determined by integrating the decay distribution obtained by using as inputs the slope 
parameters determined from the different methods of fitting, and fixing the normalisation by comparison with 
the total decay rate T'^ defined in Eq. (IIII.2I) . and evaluated with the exact form factor The 

corresponding numerical value is given in Eq. (IIII.12I) . We denote by r the ratio between the value of fso 
determined this way from the branching ratio, and the true value computed from the exact form factor, i.e. the 
reference value fso = 1.381. 
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A1 

A2 

Reference 

xVNdof 

(9.1 ±4.3)/9 

(6.9±3.6)/7 


fso/F^iO) = fs/F^{0) 

1.38 ±0.01 

1.38 ±0.01 

(1.381,1.395) 

f'so/fsO 

0.18 ±0.03 

0.23 ±0.03 

0.191 

fso/fso 

-0.03 ±0.03 

-0.33 ± 0.44 

-0.059 

fs/fs 

0 

1.91 ±6.31 

0.199 

fs/fs 

0 

-1.01 ±6.38 

-0.032 

r 

1.00 ±0.01 

0.98 ±0.04 

1 


B1 

B2 

Reference 

xVN^of 

(8.0±4.0)/8 

(6.0±3.2)/6 


fso/F^{0) 

1.38 ±0.01 

1.38 ±0.01 

1.381 

f'so/fsO 

0.18 ±0.04 

0.21 ±0.12 

0.191 

fso/fsO 

-0.03 ± 0.04 

-0.27 ±0.39 

-0.059 

fsj fsO 

1.33 ±0.40 

0.97 ±0.45 

1.010 

fs/fs 

0 

0.43 ± 12.8 

0.199 

fs/fs 

0 

0.95 ± 14.6 

-0.032 

r 

1.00 ±0.01 

0.97 ±0.05 

1 


Cl 

C2 

Reference 

xVNdof 

(3.3 ±0.1) • lOVlS 

(15.0±5.6)/16 


fso/F^{0) 

1.38 ±0.01 

1.38 ±0.01 

1.381 

f'so/fsO 

0.17 ±0.03 

0.19 ±0.03 

0.191 

fso/fsO 

-0.02 ±0.03 

-0.04 ±0.04 

-0.059 

fs/fsO 

1.07 ±0.01 

1.01 ±0.01 

1.010 

fs/fs 

0 

0.19 ±0.01 

0.199 

fs/fs 

0 

-0.03 ±0.01 

-0.032 

fs'/fs 

0 

0.01 ±0.01 

0.012 

r 

1.00 ±0.01 

1.00 ±0.02 

1 


TABLE I: Results of the different models (1-2) and methods (A-B-C), compared to the reference values, obtained from 
the Taylor expansion of the exact form factors (IIII.4I) , (1111.511 , and (IIII.BIl with the values of the parameters discussed in 
the text. The value of the ratio r does not result from the fit, but is obtained once the fit has been performed, see the 
last paragraph of Section Till Cl The column label A1 (A2) refers to the fit method A using Model 1 (2), and so on. 


D. Discussion of the results 

In order to obtain an estimate of the uncertainty attached to the coefficients of models 1-2 using methods 
A and B, we will perform fits of the models on a series of 10000 pseudo-experiments, generated by assuming 
that the data are random variables with a mean given by our theoretical model for the neutral scalar form 
factor and a standard deviation given by the relative uncertainty of the corresponding form factor measured in 
Ke 4 decays by the NA48 experiment [H, S|. We will then determine the mean and the variance of the resulting 
distribution for each coefficient of the parameterisation considered. The results are gathered in Table. ID The 
column labeled “reference” provides a comparison with the coefficients obtained from the Taylor expansions of 
the form factors given in Eqs. (IIII.5I) and (IIII.7I) . as described after Eq. (IIII.15|) . 

As shown by the comparison between Models 1 and 2 for methods A and B, the higher powers of are 
only weakly constrained. Model 1 is very rough and provides a very poor description of the charged form factor 
(modelling it as a simple constant), which explains the very bad Xmin method Cl. Only fso and /'o//so 
can be determined with a good accuracy, but there is no significant bias introduced by the fitting procedure 
with respect to the reference values. Despite of its shortcomings, method A gives good results for the neutral 
parameters. As expected, compared to method B, method C provides a much better accuracy on the neutral 
parameters since the charged ones are constrained in this method. As shown by the ratio r, both methods yield 
accurate values of fso (at the few percent level) obtained by integrating over the phase space to consider the 
branching ratio, even methods that do not attempt at describing the < 0 region correctly. This can be easily 
understood: both methods are constrained to describe correctly \Fg (s)p for small > 0 (as can be seen by 
their agreement concerning fso and /^q): but they may differ for q^ < 0 (which exhibit larger uncertainties). 
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D2 

D2 with cr[A'°4]/10 

Reference (sO, s) 

/Ndof 

(4.8±3.1)/5 

(5.6±3.3)/7 


fso/F^O) = fs/F^iO) 

1.38 ±0.01 

1.38 ±0.01 

(1.381,1.395) 

f'so/fsO 

0.20 ±0.12 

0.20 ±0.01 

0.191 

fsJho 

0.21 ±0.88 

-0.05 ±0.06 

-0.059 

f's/fs 

-14 ±234 

0.62 ± 1.86 

0.199 

fs/fs 

32 ±516 

0.00 ± 3.95 

-0.032 

0.0 - al 

0.21 ±0.16 

0.25 ± 0.03 

0.265 

fflo 

-0.83 ±2.87 

0.00 ± 0.43 

-0.045 

r 

0.92 ±0.15 

0.99 ± 0.02 

1 


E2 

E2 with cr[A'^“4°]/10 

Reference 

xVNdof 

(13.2±5.2)/14 

(12.9±5.1)/14 


fso/F^{0) 

1.38 ±0.01 

1.38 ±0.01 

1.381 

fso/fsO 

0.22 ±0.05 

0.20 ±0.01 

0.191 

fso/fso 

-0.09 ±0.08 

-0.05 ±0.02 

-0.059 

fs/ fsO 

1.01 ±0.01 

1.01 ±0.01 

1.010 

f's/fs 

0.19 ±0.01 

0.19 ±0.01 

0.199 

f's'/fs 

-0.03 ±0.01 

-0.03 ±0.01 

-0.032 

f's"/fs 

0.01 ±0.01 

0.01 ±0.01 

0.012 

flo - “o 

0.25 ±0.10 

0.26 ± 0.03 

0.265 

2 

aQ 

0.05 ±0.16 

0.07 ±0.53 

-0.045 

r 

1.00 ±0.01 

1.00 ±0.01 

1 


TABLE II: Results for methods D2 and E2, where tttt scattering lengths are htted in addition to data on the neutral 
scalar form factor (for method D) or for both neutral and charged scalar form factors (method E). The third column 
corresponds to the case where the uncertainties for the neutral scalar form factor are divided by 10. 


However, this region is very narrow (4M^o < s < 4M|) and its contribution is further suppressed by phase 
space. Therefore, the impact of this region on the estimation of the branching ratio is very small, and the latter 
is completely dominated by the region s > 4M^ where all parameterisations agree (the uncertainties reflecting 
mainly the uncertainties of the inputs and the lack of data at large g^). 

From this discussion, one thus expects that using the fit function given by Eq. (9.1) of Ref. and 
described at the beginning of Sec. |lTl will lead to similar results for the ratio r, despite the fact that this 
model complies with the expected structure of the cusp only if one imposes strong assumptions, and should 
be considered as a mere phenomenological parametrisation to reproduce data smoothly. One indeed obtains 
r = 1 ± 0.01 and x^/^dof = (8.0 ± 4.1)/9 for method A, and r = 1 ± 0.01 and x^/^dof = (8.1 ± 4.1)/8 for 
method B, illustrating once more that a smooth parametrisation of the curve above the cusp in good agreement 
with the data is enough to obtain an accurate and unbiased value for the normalisation fso- 

The outcome of this discussion is that the measurement of \Fg p allows for an accurate determination of 
fso (at the percent level), in the current experimental setting. As shown by the ratio r, the value of fso obtained 
from the computation of the branching ratio is equal (within uncertainties) to its true value for all methods and 
parameterisations considered here. Even though one has to keep in mind that this observation is done using the 
pion scalar form factors rather than the actual Kg 4 form factors, it nevertheless suggests that the fit procedure 
adopted in Ref. Q does not bias the determination of fso, and thus cannot explain the surprisingly higher value 
of fs extracted by the NA48/2 collaboration from the AT™ channel, as compared to the value for fs determined 
from the channel. 


E. Constraining the scattering lengths 

The presence of a cusp similar to the one observed in the three-body —>• 7r+7r°7r° decay suggests that it 

should, in principle, be possible to extract information on the scattering lengths from an accurate measurement 
of the differential decay rate. At leading order, the cusp is related to the difference of scattering lengths 
Qq — Qq. Going to higher orders in (i.e. including (/?g) will also add a (weaker) dependence on a§. The 
scattering lengths can be determined only once the relative normalisation of form factors involved in Fq and 
Fi is fixed, which requires the determination of the charged parameters in some way. We define two methods 
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for this purpose. Method D is exactly as method A, without including any experimental information on Oq and 
ttg in the [i- 6 - removing them from the described in eq. (IIII.16I) ]. and similarly for method E with 

respect to method C. We proceed as before, but now also htting the scattering lengths. 

We consider model 2, as model 1 yielded poor results in the previous section for method C. In order to 
discuss the potential impact future experimental improvements could have, we consider also a situation where 
all statistical errors are reduced by 10 (but the number of bins is unchanged) for the neutral channel, keeping 
the uncertainties unchanged for the charged channel. 

The results gathered in Table HIl show that the current statistical uncertainties yield a relative uncertainty 
on Oq — Oq of around 80% for D2 and 40% for E2. For D2, the charged parameters are only very poorly 
constrained, but this does not prevent the fit to be reasonable. Reducing the statistical uncertainties by 10 (for 
the neutral part) yields a significant reduction in the uncertainties, leading to a relative uncertainty on Ug — a§ 
of 27% for D2 and 10% for E2. At this level of accuracy, there is no significant bias in the value of ajj — a§ 
extracted through these various approaches. As expected, no relevant information can be obtained on Ug, due 
to the very small sensitivity of the neutral-pion channel to this quantity. To illustrate this point, if instead we 
fix the value of Og to its central value in Eq. (IIII.16|) , our results concerning the uncertainty on Og — Og and the 
quality of the fit remain unchanged. 

From this discussion, one can hope to get some information on Og — Og using model 2, should a larger data 
set become available for AT™ in the future. One has however to keep in mind that we have assumed the equality 
between the charged and neutral normalisations in the polynomials for the neutral scalar form factors in the 
case of method D, as well as the equality between the charged parameters in the polynomials for the charged 
and neutral scalar form factors in the case of method E. These assumptions are certainly reasonable considering 
the current uncertainties involved, but one might need to reassess them in the presence of more accurate data. 
In this context, it is also interesting to notice that the current result from the DIRAC experiment Ref. [2^ is 
l®o ~ ®ol = 0.253 ± 0.Oil, i.e. a 4.3% uncertainty, so that a substantial increase of the statistical sample of AT™ 
decays is needed in order to reach a comparable accuracy. 


IV. RADIATIVE CORRECTIONS TO THE K°2 DECAY RATE 

In this Section we now discuss radiative corrections, which were addressed differently in the analyses of 
the Ar^“ and A'™ channels so far. In the latter case, no radiative corrections were applied to the decay rate 
measured in Ref. This accounts for the unspecified factor Sem in Eq- dnj. It is thus natural to ask how 
much the observed 6.5% discrepancy [see Eq. (na)] in the normalisation of the form factor measured in the 
two channels is due to this correction factor. Our aim here is not to provide a complete discussion of radiative 
corrections in the AA 4 decay channels at a level of sophistication that would match the treatment of isospin 
breaking due to the difference between masses of the charged and neutral pions. We rather want to work out 
these corrections in a somewhat simpler framework, trying to reproduce a treatment of radiative corrections in 
the neutral channel similar to the one that was applied in the charged channel, in order to make the comparison 
as meaningful as possible. 


A. Treatment of radiative corrections in data 

Let us recall how radiative corrections are treated in the charged channel Si- First, virtual photon 
exchange between all possible pairs of charged external lines are considered, and the corresponding Sommerfeld- 
Gamow-Sakharov factors are applied. The corrections induced by emission of real photons are treated with 
PHOTOS [ 29143 ^ . The latter also implements the wave-function renormalisation on the external charged legs. 
The couplings of photons to mesons are treated as point-like interactions, given by scalar QED. The result is 
then free from infrared singularities. Furthermore, one neglects the contributions that vanish when the electron 
mass goes to zero, which is a sensible limit to consider for the Ke 4 decay channels. 

Apart from the Soinmerfeld-Gamow-Sakharov factors, some contributions that would arise within a more 
systematic approach, provided by the effective low-energy theory of QGD and QED for light quarks and leptons 
(34 - 1^ . as applied in Refs. [3-[ll| to the channel with two charged pions, are not considered. These include, for 
instance, all structure-dependent corrections, where the photon is emitted from the tree-level vertices or 
from internal charged lines. The outcome of such a truncated calculation is affected by an ultraviolet divergence, 
which is removed by renormalizing the coupling |I4spGp [note that in Eq. (12) of Ref. the factor (a/rr) 
has been inadvertently omitted]. 



9 a A^ 

1-m — 

4 TT 




(IV.17) 
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This same correction factor also appears in Ref. [2^, with the ultraviolet cut-off A taken equal to Mw- From 
this last reference, we also see that the factor 9/4 decomposes as 9/4 = 3 x (1/2) — 1/4 -|- 1/2 -|- 1/2, where 
the first contribution comes from the wave-function renormalisation of the three charged mesons, the second 
from the (charged) lepton wave-function renormalisation, and the last two ones from the virtual photon loops 
between the charged kaon and the charged lepton on the one hand, and between the two charged pions on the 
other hand [the remaining divergencent contributions of this type, i.e. a photon line connecting the external 
kaon to the charged-pion lines, or the charged lepton with each of the two pions, cancel pairwise]. 

In the case of the channel, we therefore expect that this factor becomes 1 x (1/2) — 1/4-1- 1/2 = 3/4. 
Since it differs from the previous one, it cannot be absorbed by the renormalisation of the same prefactor 
iKispGp as before. It seems more natural instead to absorb these ultraviolet divergences into the normalisations 
of the form factors 


rbare 


9 a 


= fs 


rh&re 

J sO 


3 a A^ 


— fsO- 


(IV.18) 


This is also more in line with the effective theory approach mentioned above, where the form factors are also 
corrected by (different) contributions from the low-energy constants Ki or Xi which are renormalised 
by the ultraviolet divergences coming from the photon loops. Using instead Eq. (IIV.17I) in both cases would 
leave a remaining cut-off dependent contribution to the AT™ amplitude. For a typical value of A = I GeV, this 
would modify Eq. (IV.52|) at the per mille level. 


B. Radiative corrections a la PHOTOS for the decay rate 


In the following, we will try to estimate the potential impact of PHOTOS on AT™ rather than pursuing 
an effective field theory approach. If we want to reproduce the analogue of the PHOTOS treatment [IH of 
radiative corrections for the AT™ decay rate, we need to consider the wave-function renormalization of the 
charged lepton and of the kaon in (scalar) QED, and the vertex correction corresponding to diagram (a) in Pig. 
m Using a Pauli-Villars regularization, and taking the photon propagator in the Eeynman gauge, we reproduce 
the expressions of Eq. (6) in Ref. [24| for the former. In order to evaluate and discuss the contribution from 
diagram (a) in Fig. [21 we chose to describe the tree-level Ar°° vertex as 

= [F™(pi + A2 )m + R°%k - Pi - P2)^] . (IV.19) 

with constant form factors [Bose symmetry forbids a contribution of the form G°°(pi — P 2 )/i with constant], 
so that the lowest-order amplitude reads 

A(a:°^) = -^^Vusuipeh^ii - j 5 )y{p^) [a°°(pi +P2)t^ + R°°{k-pi -P2)/.] • (iv.20) 

In the limit rUe —>■ 0, we obtain 


A{K°Jma)) = e^Ao{K°J) 




{{k - pef;MK,me) - A{k ■ pe)Ci2 ((fc - Pe)^; me)] 


R°°{P1+P2)u X u(pe)7‘'(l - 75 )v(Pi,J 


Mk 

1 / 1 

3^ (^31n^ -b - ) -4(fc-pe)Gii {{k- Pef-,MK,rne) +2J ({k - Pef; Mx^me) 


(IV.21) 


The various loop functions occurring in this expression are defined in Appendix|5] Adding to it the wave-function 
renormalizations on the charged external lines gives the following result, in the framework adopted here, for 
the radiatively corrected amplitude [rrij denotes a small photon mass, introduced as an infrared regulator, to 
be sent to zero once an infrared-safe observable has been constructed]: 
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(47r)2 


3, A2 
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-4(fc.pe)Ci2((fc 


PeY\MK,m^) 


X 


+P2)u X u(pe)7‘'(l - 75)v(p^J 

-4(fc • Pe)Cii {{k - Pef] Mk, TOe)] . 


(IV.22) 


We make a few comments about this result: 


• Although the result pV.21|) holds in the Feynman gauge ^ = 1, we have also computed the wave-function 
renormalisations and A( K^^ma)) in an arbitrary linear and covariant ^-gauge, and we have checked that 
the final result pV.22l) does actually not depend on the gauge-hxing parameter 

• In order to reproduce the analogue of the PHOTOS treatment [l^, of radiative corrections for the AT™ 
decay rate, one needs to add the emission of soft photons from the charged external lines, diagrams (a) 
and (c) of Fig. [3] [in the Kjl case, there are two more diagrams where the photon is emitted from the 
charged pion lines] so that the result is free of infrared singularities at order a. These corrections will be 
discussed later on. At this stage, we simply note that the infrared-divergence of Eq. (IIV.22|) is equal to 

X Inm.^ [-4 - 2(k ■ pe)T(k,pe)] , (IV.23) 

with the function r(pi,p 2 ) defined in Eq. (|B.20p . 

• The factor | f In discussed in Eq. (IIV.18|) is indeed to be found in Eq. (IIV.22I) , provided that one sets 

to zero. The only remaining contribution comes from which is proportional to /so at this level. 
This indeed corresponds to the situation considered in Ref. [^. In the absence of radiative corrections, 
the form factor [or i?'* in the charged channel] does not contribute to the decay distribution for 
me = 0. In this case, one may as well take RP^ = 0 from the beginning. But once radiative corrections 
are switched on, taking me = 0 or RP^ = 0 are no longer equivalent options. As shown by the second 
contribution in Eq. (Iiy.22[) . there is a correction to that is induced by i?™, and this contribution is 
not considered in Ref. |24|. and is hence also missing in Ref. (3^ . 

• At lowest order, and in the isospin limit, one has [s^ 


i?o° = 


pOO 

1 1 

pOO 

'2 



2 

[ Se-Ml\ 


3 

2 

/3 Se-M],)\ 


(IV.24) 


Actually, as shown in the second expression, the vertex in the diagram (a) in Fig. [2] only accounts for the 
contribution = (2/3) • The second factor comes from the diagram (/) in Fig. [2| 


C. Additional non-factorizable radiative corrections to the A /4 decay rate 

After these preliminary remarks concerning the PHOTOS-type treatment of radiative corrections in the 
and A /4 channels, let us now address radiative corrections in the channel with two neutral pions in a 
somewhat more systematic manner. This will allow us to estimate the size of the radiative corrections that are 
not included in the experimental analysis, as described in the previous section. We keep on considering the 
limit where me vanishes, so that in the absence of radiative corrections the amplitude reads simply 

■P-oiKei) = i-^Vusn{pe) {^1+ h) (1 - X + 0{me). (IV.25) 

For our purpose, it is convenient to distinguish between two types of radiative corrections, that we call factoris- 
able and non-factorisable. Factorisable radiative corrections are defined by the contributions where both ends 
of the virtual photon line connects to a charged mesonic line or to the vertex with the leptonic current, or when 
both ends connect to the electron line. These factorisable contributions will not modify the structure of the 
matrix element, but will change the form factors and i?°°. We find it convenient to express them as [dll] 

Q pOO 

Afact(A//) = i—^Vusn{pe) {^1+ ^ 2 ) (1 - 75)v(p,.) X X ^/^^/^ + 0{me), 

V2 


(IV.26) 
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where we have factored out the wave-function renormalization factors computed in QED for Zg, and in scalar 
QED for Zk- 


= 1 
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In contrast to the preceding subsection, we use now dimensional regularization, with the minimal subtraction 
of the combination 


A = 


1 




^-i(ln(47r)+r'(l) + l) 


(IV.28) 


It is understood that in Eq. (|IV.26() now includes all the remaining factorizable photonic corrections, 
together with the contributions from the low-energy constants Li [s^ and Ki [s^. These will take care of the 
UV divergences due to the meson loops and to the photon loops, respectively, so that the product is 

actually UV finite. It however inherits the infrared divergence contained in yjZ k- 

Let us next consider the non-factorisable contributions. As far as the corrections to are concerned, 
they are represented by the diagrams shown in Eig. [21 On finds that the contributions coming from the diagrams 
(6) and (c) are proportional to the lepton mass, and thus vanish in the limit mg —^ 0. There are therefore only 
three diagrams to compute in this approximation. Consistently dropping terms that vanish as mg —>■ 0, one 
finds that 


A(Klim.a)) 


Ao{Kll) X [4(fc • pg)C ((fc- 
- J {{k - pef']MK, mg) - 4A 


Pef] Mk, me) - 2{k ■ Pe)Ci2 {{k - PeY', Mk, me) 
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A[K^elme)) = Ao{K°e!) x y 


J {{k - Pi - P2Y,m^, Mk) - 2A 


1 

IOtt^ 



and 


A(Kiim)) 


Aq{KII) X y [2{PI+P2fc ((pi -bp2)^, {k-pi-p2f]MK,MK) 

-(pi +P2)^C'ii ((pi -bP2)^, {k-pi -P2)^; 
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(IV.31) 


Apart from the change of regularization, the expression for Ah fSIaU in Eq. (IIV.29P reproduces the one of 
Eq. (|IV.21|) obtained previously, provided one takes = (2/3) • A™, as discussed at the end of Sec. IIVBI 
and makes use of the identities 

{k-pe)Cxi{{k-pef-,MK,me) = --b ^ J ((fc - pg)^; M/f, mg)-f ..., 

1 / M^ A 1- 

MlCii{(k-pef-,MK,me) = f l-fln-^j -b-J ((fc-pg)^;Mif,mg) 

-{k ■ Pe)Ci2 ((fc - pg)^; Mk, me) + ..., (IV.32) 


where the ellipses denote terms that vanish in the limit mg —>■ 0. 

Adding up the contributions discussed so far, one obtains an expression for the radiative corrections at 
order 0{a) that still contains both infrared and ultraviolet divergences. The latter will be taken care of by the 
contributions from the counterterms Xi introduced in Their contribution reads 


Al(ifO°;cts) = A(ihgOO) X - ive) . 


(IV.33) 


The low-energy constant Xi is not renormalized, whereas Xq = Ag(p) — 5A. Collecting the divergent pieces 
from the various contributions leads to [we recall that at this stage F^^^Zk has already been made UV-finite 
through the contributions of the low-energy constants Li and Ki] 


A(Ag° 4 °; UV-div) 


e^Ao{K°J) X 
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FIG. 2: The set of one-loop virtual photon exchange diagrams of the non-factorisable type to consider for the K°2 decay. 


which vanishes as it should. 

As to the infrared divergences, collecting the IR-divergent pieces contained in the contributions computed 
so far, one obtains 


A{K°2; IR-div) 



Ao{K°°) X In TO-), 


+ i-‘^){k-Pe)T{k,Pe) , 
y/Z^ yf Zk Ha) 


(IV.35) 


with the function r(pi,p 2 ) defined in Eq. (IB.201) . Besides the wave-function renormalisations, such divergences 
only arise from the contribution of C ({k — pe)^; MK,rne) in A{K^2i Notice that this infrared divergence 

coincides with the one of the result (I1V.22I) . given in Eq. (IIV.23I) . The construction of an infrared-safe observable 
at order 0(a) requires also to consider the process with the emission of one soft photon. The corresponding 
differential decay rate is given by 


^ X (27r)'‘<5‘‘(fc-pi-p2-p^-Pa-9) 
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FIG. 3: The one-photon emission contributions to the decay. 


The amplitude for the radiative decay 7 r° 7 r°e^ I'e 7 can be expanded in powers of the photon 

energy, 




(IV.37) 


The Low approximation consists in keeping A-i(K^ 4 j) alone. This is enough in order to study the emission of 
only soft photons and to discuss the issue of infrared divergences. Explicitly, one has [q^ is the momentum of 
the emitted (real) photon, e*{q) the corresponding polarization vector] 


= eAoiK^,2) 


00 ^ I PC ^*( 9 ) k-£*{q) 


\Pf q + ^ k-q- 
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Then 
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One may then perform the integration over the undetected soft photon. In the soft-photon approximation, the 
photon momentum in the delta-function of the phase-space integration is neglected, and one takes 


^rsoft(^oo^) ^ {2TT)^S'^{k -Pi-P2-Pi- Pv) X 


d^pi (fp2 cPpi (fpi, 


" 2 (2i2|,1 > 


(27r)32£;i (27r)32E;2 {2TTf2Ei {2Trf2\p^ 
2 


= dro{K°!) X (-e^) / 


d^q 


-2 


\q\<AE (27r)32|^ 
k-pt 


Ml 






(IV.40) 


Expressions for the corresponding integrals can be found in ip. The integral is limited to photon energies 
1^ below the experimental detection threshold AE in the kaon rest-frame. As far as the infrared divergences 
are concerned, one has 


X ^ X Inm-, [2 + 2 + 2(fc • pMk.Pe)] 


(IV.41) 
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Therefore, the contributions proportional to Inm^ cancel in the sum dr{K^2) + dT{K^2-y)- For later convenience, 
we rewrite Eq. (IIV.40I) in a way that explicitly displays the IR-singular part: 

= dVoiK^J) X ^ [2 + (fc . PeMk,p,)] X 2 In ^ (IV.42) 

We then add the contribution 

9 

777 

MK°2) X [2+ (A:.pe)r(fc,pe)] x 21n^ (IV.43) 

to the amplitudes involving virtual photons, such as to make them infrared finite. To this end, we define the 
function 

C {{k-pef;MK,me]AE) = C {{k - pef', MK,me) + 3 ^ ^ (IV.44) 


D. Discussion on radiative corrections for K^2 


We can now add the virtual and real contributions described up to now which should be involved in a 
PHOTOS-like treatment of this decay. We include them as a correction of the form (1+5em) to the determination 
of the form factor from the measurement of the branching ratio. To this end, we compute the total decay rate 
including the soft photon emission 

ptot = r{K°!l) + (iv.45) 

where includes corrections at first order in the fine-structure constant a, and write it in terms of the 

decay rate ro(i^g 4 ) without radiative corrections in the form 

= ro(if° 4 °) X (1 + 25em) ■ (IV.46) 


Let us first discuss the corrections computed in Subsection B above. In order to obtain a result that is as 
close as possible to the treatment of radiative corrections in the channel, we absorb the UV-divergent factor 
of Eq. (IIV.22P in and take i?™ equal to zero. Then r{K^2) is computed by performing the phase-space 
integration of 

A(iF° 4 °) [A iK°!) + 2e^AA{K°J)] , (IV.47) 


where 
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We take AE = 11.7 MeV, the value corresponding to the NA48/2 experiment Q, for the real-photon detection 
threshold in the kaon rest frame. This gives then 


Sem = 0.018. 


(IV.49) 


This value has the expected size. Moreover, it goes into the right direction, in the sense that it reduces the 
discrepancy in Eq. (II. 2^ from 6.5% to 4.6%. 

As a test of the stability of the result (IIV.49|) we may also evaluate the non-factorizable radiative cor¬ 
rections corresponding to all the diagrams in Fig. [2] This amounts to taking the expressions in Eqs. (IIV.29L 
(IIV.30() . and (IIV.31I) for the evaluation of r(A"g 4 ) (let us stress again that these equations have been obtained 
in a regularisation scheme differing from the one discussed in Sec. lIVBl) . We absorb the ultraviolet divergences, 
as well as the contribution (IIV.33|) into in order to build a quantity both UV and IR finite. Constant 
terms have been discarded as they could be included in the contribution of the counterterms Ai and Xq. The 
resulting expression for AA{K°2) then reads 
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+ {pi + P2fC {{pi + P2f, {k-pi- P2f] Mk, Mk) 

-^(pi +P2fCii ((pi +P2f, (k-pi - P2 )^-,Mk,Mk) 
{{k-P i-P 2f\rn^,MK) , 

instead of the expression in Eq. (IIV.48I) . For AE = 11.7 MeV, we obtain now 


(1V.50) 


5em = 0.017. (1V.51) 

This value is quite close to the one obtained in Eq. (I1V.49I) . so that in the present case the treatment of radiative 
corrections a la PHOTOS seems to yield stable results even after the inclusion of non-factorizable contributions. 


V. SUMMARY AND CONCLUSION 

The present study is devoted to isospin-breaking effects in the semileptonic decay of the charged kaon 
into two neutral pions, —>■ Ve ■ Because of the smallness of the electron mass and of the limited 

experimental precision, this decay can be described in terms of a single form factor. This form factor also 
occurs in the description of the decay into two charged pions, > 7r“'"7r“e^ r'i, and up to isospin-breaking 

contributions, the two determinations should agree. The present study focuses mainly on two aspects related to 
this issue: i) to ascertain quantitatively to which extend the phenomenological parameterizations used in order 
to analyse the data could impinge on the resulting value of the normalization or on the shape of the form factor 

measured in the decay —>■ Tr^ir^e^ r'J, and ii) to obtain a quantitative estimate of the radiative corrections 
to the total decay rate, which again might affect the normalization of the form factor. 

Concerning the first issue, we have considered the form factors of the pion as a case study. As a first 
step, we have discussed the structure of the form factors, and their properties linked to the presence of a 
cusp, using exact expressions of the form factors valid up to two loops in the low-energy expansion. We have 
clearly established that the phenomenological parameterizations used in order to analyse the data did not 
agree with the general properties that can be infered form these exact expressions. In a second step, we have 
generated pseudo-data from these form factors, that we have then analysed using several phenomenological 
parameterizations. The outcome of this study is that the determination of the normalization of the form factor 
is actually not sensitive to the parameterizations used. As a side product, we see that the higher orders in 
the Taylor expansion of form factors are not accurately determined by a direct fit to simplified (polynomial) 
formulae. Although our study was carried out for the scalar form factor of the neutral pion, we expect that 
the conclusion also holds for the AT™ form factor. We have also considered the possibility to constrain the tttt 
S'- wave scattering lengths from the measurement of the decay distribution. We have found that, unfortunately, 
with the sample of events presently available, the statistical uncertainties remain large. A statistical sample 
comparable to the one available in the channel would be required in order to reach a precision close to 

that obtained by the Dirac experiment. 

The second issue addressed in this paper consists in radiative corrections. We have determined the cor¬ 
rection factor Sem to the total decay rate in Eq. (11.21) . In order to make a meaningful comparison with the value 
of the normalization of the form factor extracted from the A'^~ channel, we have used a simplified framework, 
including only those corrections that were also included in the latter case (one-loop photonic corrections on the 
wave functions and tree-level vertex). Our result Sem = 0.018 leads to the replacement of Eq. (11.21) by 

\ = 1.046(10). (V.52) 

fslKir] 

Note that the error bar in this equation is purely from experimental origin, and does not include the systematic 
uncertainties from the methods used for the evaluation of radiative corrections in both channels. Such additional 
uncertainties can stem, for instance, from the regularisation dependence of the PHOTOS(-like) treatment of 
radiative corrections, and from neglecting the dependence in the cut off A discussed in Sec. IIV Al 

We have also considered additional photonic corrections estimated within a different regularisation 
scheme, and we have found that they do not modify the previous estimate in a significant way. A few comments 
are in order: 

• The analysis of radiative corrections we have performed provides an adequate estimate of the global factor 
Sem that modifies the total decay rate. It need not be suitable for an analysis of radiative corrections to 
the phase-space distribution itself. 
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• Other isospin-breaking corrections, among them factorizable exchanges of virtual photons, but also effects 
due to m„ — or to the mass differences between charged and neutral pions and/or kaons, are not 
covered by our analysis. They could affect the normalization of the form factors measured in the two 
channels in different ways. A mode elaborate study is needed in order to reach a quantitatively meaningful 
interpretation of the result in Eq. (IV.521) . 


• At lowest order in the chiral expansion, these additional isospin-breaking corrections are given by Eq. 
m- For R = 35.8(1.9)(1.8) [3, and adding errors in quadrature, we obtain 






= 1.042(3). 


(V.53) 


LO 


In view of the value given in Eq. (IV.52|) , the corrections from higher orders to this ratio should therefore 
be small. 


• Conversely, using the relation m in regard to the result (IV.521) . and discarding yet to be computed 
corrections to the former, we obtain R = 32/^g. 

The discussion of the radiative corrections presented here is clearly only a first step. In view of the 
statistical accuracy of the data, a full model-independent calculation of these corrections in the neutral as well 
as in the charged channels is certainly mandatory before a definite conclusion can be reached concerning the 
observed difference in the normalisation of the form factors between the neutral and the charged channels. This 
task is clearly beyond the scope of the present note, and is left for future work. 
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Appendix A: Properties of the functions A“(s) 


In this Appendix, we wish to summarise the properties of the functions A'“(s) that are needed in the 
discussion of the cusp in Section im Let us start with the functions A'°(s) and K^{s), defined by dispersive 
integral as in Eq. (III.261) . with sthr = 4M^o and 371 


fco(s) = 


= ^Lois). 


k^,{s) = 

kJis) = 




1 

87r 
1 

87r cr(s -I- 4A,r) 


1 - 


CTo(s) 


Lois), fc°(s) = 


Mi 


IGtt scto(s) 


Llis), 


Lis + AA^), fcj(s) = 


Mi 


IGtt stTo(s) 


-|- dAjr) 


(A.l) 


The definitions of the functions cto(s) and Lo{s) for s > AM^o can be found in Eqs. (III.8|) and (III.11|) . respectively. 
The function cr(s) is also to be found in Eq. (III.81) . whereas L{s) is defined as 


[ ln(^)^L(s) [s>AMi] 

L{s) = I , (A.2) 

[ln(^)^T + *7r [4M2„<s<4M2] 


according to the definitions (III.141) and (III.15|) . For s > AM'^q, the functions A:°(s) and fc^(s) are real and 
smooth. In the same range of s, the functions K^is) and K^is) have smooth real and imaginary parts, with 
ImA'°(s) = fc°(s), ImA',^(s) = fc^(s). Finally, the functions K°is) for n > 1 can be expressed in terms of 
jo(s) = A'g(s). The explicit expressions and their derivation were given in 

There is not much to add as far as the functions AV(s) are concerned: it is sufficient to replace everywhere 
in Eq. (IA.1|) M^o by the charged pion mass M^, and hence (To(s) by cr(s), and Lo{s) by L{s). In the dispersive 
representation (III.261) . the integration starts at Sthr = AM^. In the case of Ko{s) = the decomposition 

(III.121) then follows from (|II.28p by noticing that fco(s)/cr(s) is a constant, and that [22, [sj 


Re J(s) 


1 

IfiTT^ 


[2-I-(t(s)L(s)] [s>0]. 


(A.3) 
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For the remaining functions Kn{s), it is most convenient to use their expressions in terms of J(s). One then 
finds 


= L{s)'^ — TT^ 


= L[s){L^{s) - 27r^) 


1671^/^2 = — 4 

Ml TT^ 


S(t(s) 2 


leTT^/irf' = -27r^^ 

(t(s) 

1677^^2^^ = —27rL(s)(7(s) 

3 s a^{s) 


(A.4) 


One may check that all these functions are real and smooth for s > 4M^o [actually, for s > 0]. 
The functions Kl^{s), n = 1,3, are defined by 

f'xt \ s f°° dx k^ix) 

Kis) = - /- 

TT 74M2 ^ X — S — lij 


(A.5) 


with [the quantities appearing in these formulae are defined in Ref. see Eqs. (2.13), (4.13), (4.16), and 
(4.17) therein] 


kti^) = 


1 1 


Stt scto(s) - 


A^/^(t_(s))£^(s) - A^/2(4 (s))/:+(s) 


^ 2 ( 5 ) = -^^is)aois)Lo{s), 

Stt 

= il;;S^ 


(A.6) 


The three functions kl^(s), n = 1, 2,3, are real and smooth for s > Ml, and they become purely imaginary for 
4A7 ^o < s < 4A£^. The functions fc^(s) = kl{s)/a{s) are then smooth in the range 4M^o ^ s < M]^, so that 
one obtains 


iF)?W(s) = - 


fcn(g) 

cr(s) ’ 


(A.7) 


and 


ReiFf,3(s) 


[s > AMI] 




Ks) 


(A.8) 


Re^f.3(s) + ^*i.3('S) ^ ^ 4M2] 


Since analytical expressions for Re A'f 3 ( 5 ) are not available, one has to use the integral representation given in 
Eq. (IA.5|) for numerical applications. 

Finally, there remains to discuss the function /C“(s) whose discontinuity along the real s axis for s > AMl 

reads 


r(s) = 


1 Ml ct(s) 1 


(ct(s) - cro(s)) A^/^(t_(s))£_(s) - (cr(s) + CTo(s)) A^/^(t+(s))£+(s)j . (A.9) 


IGtT A^n- l7o(s) S 

The function k^{s)/a{s) is real and smooth for s > 4M^o. Thus one infers 

k^{s) 


/C“W(s) = 


a{s) 


(A.IO) 


and 


£"M(s) = < 


ReK.^{s) 


[s > AMI] 


t{s) 


(A.ll) 


Re/C"^(s) + ^A:“(s) [4MA <s< AMl] 
cr(s) 
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Appendix B: Loop functions 


The computation, in Section 113 of the diagrams describing the virtual photon corrections involves a 
certain number of loop functions, that are briefly discussed here, in order to make the calculation in section Hvl 
self-contained. 

At the level of the two-point one-loop diagrams, one has 


J (p^;mi,m2) 


1 

(27r)4 {P — ml)[{i — pY — m^] ’ 


(B.l) 


and 


Jf! [p ;mi,ni2) = ~ J 


1 f 


{2 ttY (P — — pY — mi] 


Et 

2p' 


2 [(p^ + ml — mi) J (p^; 7711 , 7712 ) -f iA{ml) — iA{ml)] , 


where A{m‘^) is related to the tadpole graph, 


A{rrY) = 

Other useful relations are 
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- 1 771^ 


_ _ 1 777? In ^ — 777? In ^ 
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2 2 
rrii — 777-2 
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J{rrY] 0,777) = —2A -|- 


IGtt^ 


1 — In 


777 


21 


(B.2) 


(B.3) 


(B.4) 


(B.5) 


Explicit expressions of the function J(p^; mi,m 2 ) can be found in Ref. [s^. Moreover, the link with the functions 
J(s) and Jo(s) encountered in Section|n]is given by J(s) = J(s; M^^., M^r) and Jo(s) = J {s', M^o, M^o). 

As far as the three-point one-loop functions are concerned, one has 
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1 r dH 
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where p\ = mf. Explicitly, one has [A(a;, y, z) = x'^ + — 2xy — 2xz — 2yz\ 
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Ci2 {{pi - P2f,pl;'rni,m2) = 
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In the case where also pi = m|, these expressions simplify further, and one obtains 
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For pI = ml, the function C (t; mi, m2) itself reads [ 1 , [T^ 


1071 ^ X{t,ml,ml) m\ — ml m\ 

1 ^ 


C {t;mi,m 2 ) = 
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The two roots of Py are then given by 


y± = a±b. 


with 
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In the case under consideration, we have p)* = k>^, P2 = Pe, mi = Mk, m2 = me, with = M]^, pi = 
t = [k — pe)^ > 0 . Then, X{{k — pe)^,M|-, ml) > 0 , and 


“ 2 


1 - 


{Mk - me){MK + me) 
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> b= Y^XP^{t,M]^,ml), 
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with a < 0 and |a| > | 6 |. Therefore, 

11 / m‘ 

C(t;mi,m2) = ^-Inl — x In 
1 1 



y+ - 1 _ y_ 
y- -1 y+ 


1287 r 2 bt 
1 1 
647 r 2 bt 
1 1 
327 r 2 bt 


[ln^(l - y+) - ln{-y+) - ln^(l - y.) + ln^(-y_)] 

ln(l - p+)ln(l - y_) - ln(-p+) ln(-p_) - 21 n( 26 )ln (— - 

V y+ 


Li2 


y+ -1 

2 b 


- Li: 


\ 2 b) 


(B.17) 







































28 


Notice that 

?/+ — 1 ?/„ — a + b p1+ p^ — t + 2tb 

y_ — 1 — b'^ — a — b pi+pi — i — ‘2‘tb ’ 

so that the infrared divergent piece of C is given by 
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Let us now consider the case where p^ m 2 , but with pi = ml as before. Going through the same steps 
as in the previous case, one obtains 
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with = yp1 + (1 — y)p 2 - As before, one has 
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The case under consideration here corresponds to p^ = k, pi = = M^, whereas P 2 = {k — pi — P 2 )^, 

pI = Se, t= {pi -P 2 Y = Stt, with ml < Se < {Mk - AMl < < [Mk - ■ Then X{t,p1,pl) = 

A(s^, Se, Mj.) > 0, so that b is real, with |a| > |6|. On the other hand, one has 
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In the case at hand, this gives a = 1/2, and b^ = (stt — 4M|.)/(4s7r) < 0, so that 
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Then one obtains 
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